COSMOS Photometric Redshifts with 30-bands for 2-deg 2 
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ABSTRACT 

We present accurate photometric redshifts in the 2-deg 2 COSMOS field. The redshifts are 
computed with 30 broad, intermediate, and narrow bands covering the UV (GALEX), Visible-NIR 
(Subaru, CFHT, UKIRT and NO AO) and mid-IR {Spitzer/IRAC). A X 2 template-fitting method 
(Le Phare) was used and calibrated with large spectroscopic samples from VLT-VIMOS and 
Keck-DEIMOS. We develop and implement a new method which accounts for the contributions 
from emission lines ([O II], H/3, Ha and Lya) to the spectral energy distributions (SEDs). The 
treatment of emission lines improves the photo-z accuracy by a factor of 2.5. Comparison of the 
derived photo-z with 4148 spectroscopic redshifts (i.e. Az — z s — z p ) indicates a dispersion of 
(Ta 2 /(i+z b ) = 0.007 at < 22.5, a factor of 2 — 6 times more accurate than earlier photo-z 
in the COSMOS, CFHTLS and COMBO-17 survey fields. At fainter magnitudes i£ B < 24 and 
z < 1.25, the accuracy is <J& z /(i+z s ) — 0.012. The deep NIR and IRAC coverage enables the 
photo-z to be extended to z ~ 2 albeit with a lower accuracy (caz/(i+2 3 ) = 0.06 at i~^ B ~ 24). 
The redshift distribution of large magnitude-selected samples is derived and the median redshift 
is found to range from z m = 0.66 at 22 < i\ B < 22.5 to z m — 1.06 at 24.5 < i\ B < 25. At 
i^ B < 26.0, the multi-wavelength COSMOS catalog includes approximately 607,617 objects. The 
COSMOS-30 photo-z enable the full exploitation of this survey for studies of galaxy and large 
scale structure evolution at high redshift. 

Subject headings: galaxies: distances and redshifts — galaxies: evolution — galaxies: formation 
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1. Introduction 

Photometric redshifts (hereafter photo-z) are 
an estimate of galaxy distances based on the ob- 
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served colors (Baum 1962). This method is ex- 
tremely efficient for assembling large redshift sam- 
ples for faint galaxies. Despite having a lower 
accuracy than spectroscopic redshifts (hereafter, 
spectro-z), photo-z have the advantage of signifi- 
cantly improved completeness down to a flux limit 
fainter than the spectroscopic limit. Deep photo-z 
samples such as COMBO-17 (Wolf et al. 2003), 
CFHTLS (Ilbcrt et al. 2006), SWIRE (Rowan- 
Robinson et al. 2008), and COSMOS (Mobasher 
et al. 2007) contain more than 1,000,000 galaxies 
and go as faint as i ~ 25 with a relatively small 
amount of telescope time. 

Typical photo-z with an accuracy of <7a z /(i+z b ) ~ 
0.02-0.04 (Az widely used to study 

the evolution of galaxy stellar masses and lumi- 
nosities (e.g. Fontanaetal. 2000; Wolf et al. 2003; 
Gabasch et al. 2004; Caputi et al. 2006; Arnouts 
et al. 2007), for angular clustering analysis (e.g. 
Heinis et al. 2007; McCracken et al. 2007), to 
study the relation between galaxy properties and 
environment (e.g. Capak et al. 2007), to trace 
large-scale structures (Mazure et al. 2007; Scoville 
et al. 2007) and to identify clusters at high red- 
shift (Wang & Steinhardt 1998). Photo-z are also 
necessary for dark energy and dark matter weak 
lensing studies to separate foreground and back- 
ground galaxies and to control systematic effects 
such as intrinsic shape alignment, shear shape 
correlation, and the effects of source clustering 
(Peacock et al. 2006). All of these applications 
require strict control of systematic effects in the 
photo-z estimate. An efficient way of identifying 
and removing systematics is to calibrate photo-z 
on a spectroscopic sample. The most common 
methods of calibration are neural network meth- 
ods (e.g. Ball et al. 2004, Colister et al. 2004), a 
calibration of the color-z relation (Brodwin et al. 
2006; Ilbert et al. 2006), or a reconstruction of the 
SED templates (Budavari et al. 2000; Feldmann 
et al. 2006). 

As is also true for spectroscopic redshift mea- 
surements, photo-z accuracy depends on spectral 
coverage and resolution. It is also degraded for 
sources with a low signal-to-noise ratio (Bolzonella 
et al. 2000). The Balmcr and Lyman breaks con- 
tain much of the photo-z information, so accu- 
racy is lower in redshift ranges where these fea- 
tures are not well sampled by the filter set. As the 
photometric accuracy degrades, it becomes more 
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difficult to constrain the positions of these fea- 
tures, leading to lower accuracy. This creates a 
dual dependency of photo-z accuracy on magni- 
tude and redshift which is defined by the survey 
design (exposure time, filter choice, and calibra- 
tion accuracy). For this reason, any photo-z sam- 
ple should be extensively tested and characterized 
in the same way a spectroscopic sample would be. 
This analysis is necessary in order to identify the 
redshift /magnitude ranges over which the photo-z 
can be trusted and used for a given scientific ap- 
plication. 

This paper presents a new version of the pho- 
tometric redshifts for the Cosmic Evolution Sur- 
vey (COSMOS: Scoville et al. 2007) and an anal- 
ysis of their accuracy. COSMOS is the largest 
Hubble Space Telescope (HST) survey ever under- 
taken - imaging an equatorial 2-deg 2 field to a 
depth of If814w = 27.8 mag (5 a, AB). The COS- 
MOS field is equatorial to ensure visibility by all 
ground-based astronomical facilities. This project 
includes extensive multi-wavelength imaging from 
X-ray to radio (XMM, Chandra, GALEX, Sub- 
aru, CFHT, UKIRT, Spitzer, VLA). New ground- 
based NIR data (McCracken et al. 2008; Capak 
et al. 2008), Spitzer-IRAC data (Sanders et al. 

2007) , and medium/narrow band data from the 
Subaru Telescope (Taniguchi et al. 2007; Sasaki 
et al. 2008; Taniguchi et al. 2008; Capak et al. 

2008) greatly improve the previous photometry 
catalogue (Capak et al. 2007). These new data 
are used for the photo-z derived here, yielding a 
factor of 3 higher accuracy than the first release 
of COSMOS photo-z (Mobasher et al. 2007). 

The COSMOS data are presented in §2. The 
technique used to estimate the photo-z is pre- 
sented in §3. In §4, we quantify the photo-z accu- 
racy as a function of the magnitude and redshift. 
In §5, we provide the photo-z distribution of the 
i + selected samples. Specialized photo-z for X- 
ray selected sources, AGN, and variable objects, 
are discussed in a companion paper (Salvato et al. 
2008). 

Throughout this paper, we use the standard 
WMAP cosmology (O m = 0.3, tt A = 0.7) with 
h = 0.7 and h = H /100 km s" 1 Mpc" 1 . Mag- 
nitudes are given in the AB system. 



2. Data 

Compared to the previous optical/near -infrared 
catalogue (Capak et al. 2007) the new photom- 
etry implements 14 new medium/narrow band 
data from the Subaru Telescope, deep ground- 
based NIR data (J and K bands) and Spitzer- 
IRAC data. The spectroscopic sample used to 
calibrate/test the photo-z is 10 times larger at 
i\ B < 22.5 than that of Mobasher et al. (2007). 
The spectroscopic sample is supplemented with 
faint infrared selected sources and a deep faint 
spectroscopic sample at z > 1.5. Hereafter, we de- 
tail the photometric and spectroscopic data used 
to measure the photo-z. 

2.1. Photometric data 

Fluxes are measured in 30 bands from data 
taken on the Subaru (4200-9000A), CFHT (3900- 
21500A), UKIRT (12500A), Spitzer (3.6-8 urn) and 
GALEX (1500-2300A) telescopes. We refer to Ca- 
pak et al. (2008) for a complete description of the 
observations, data reduction and the photometry 
catalogue. The equivalent width and the effective 
wavelength of each filter are listed in Table [T] The 
sensitivities are given in Capak et al. (2008) and 
in Table 1 of Salvato et al. (2008). Table Q] indi- 
cates the fraction of sources detected in each band 
with an error less than 0.2 mag for a selection at 
i + < 24.5, i + < 25 and i + < 25.5. We summarize 
below the datasets used in this paper. 

Ultraviolet: Very deep u* band data were 
obtained at the 3.6m Canada-France Hawaii 
Telescope (CFHT) using the Megacam camera 
(Boulade et al. 2003). The u* band data were 
processed at the TERAPIX data reduction cen- 
teJE The u* band data cover the entire COSMOS 
field and reach a depth of u* ~ 26.5 mag for a 
point source detected at 5cr. The u* band im- 
ages are also used as priors in the measurement of 
FUV (1500A) and NUV (2300A) fluxes in order 
to ensure a proper deblending of sources in the 
GALEX images (Zamojski et al. 2007). GALEX 
fluxes are then extracted using the EM-algorithm 
(Guillaume et al. 2006). They reach a depth of 
FUV ~ 26 mag and NUV ~ 25.7 mag. 

Optical: The COSMOS-20 survey (Taniguchi 
et al. 2008) entailed 30 nights of observation at 
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the Subaru 8.2m telescope using the Suprime-Cam 
instrument. The observations are complete in 20 
bands: 6 broad bands (Bj, Vj, g + , r + , i + , z + ), 
12 medium bands (7A427, 7A464, 7A484, 7A505, 
7A527, 7A574, 7A624, 7A679, 7A709, 7A738, 
7A767, 7A827) and 2 narrow bands (JVB711, 
NB816). 

Near- infrared: The catalogue includes deep J 
and K band data obtained using the WFCAM and 
WIRCAM wide-field infrared cameras on UKIRT 
and CFHT, respectively. The NIR data reduction 
is detailed in Capak et al. (2008) and McCracken 
ct al. (2008). The data reach J ~ 23.7 mag and 
K <~ 23.7 mag for a point source detected at 5a. 

Mid-infrared: Deep IRAC data were taken 
during the Spitzer cycle 2 S-COSMOS survey 
(Sanders et al. 2007). A total of 166 hr were dedi- 
cated to cover the full 2-deg 2 with the IRAC cam- 
era in 4 bands: 3.6/im, 4.5/jm, 5.6/im and 8.0/im. 
Source detection is based on the 3.6/7777 image and 
the fluxes were measured in the four IRAC bands 
using the "dual mode" configuration of SExtrac- 
tor. The IRAC catalogue is 50% complete at l/i Jy 
at 3. 6/7777 (771,3. e Mm ~ 23.9 mag). 

All of the imaging data were combined to gen- 
erate a master photometry catalogue (Capak et al. 
2008). Photometry was done using SExtractor in 
dual mode (Bertin & Arnouts 1996). Source de- 
tection was run on the deepest image (z + ~ 26.2 
for a point source detected at 5a). For the UV- 
NIR data, the Point Spread Function (PSF) varies 
from 0.5" to 1.5" from the K to the u* images. 
In order to obtain accurate colors, all the images 
were degraded to the same PSF of 1.5" following 
the method described in Capak et al. (2007). The 
final photometry catalogue contains PSF matched 
photometry for all the bands from the u* to the K 
band, measured over an aperture of 3" diameter 
at the position of the i + band detection. For the 
FUV and NUV data, we transformed the total flux 
provided for the GALEX-FUV and -NUV coun- 
terpart by multiplying it by a factor 0.759. This 
factor is the fraction of the flux observed in opti- 
cal into a 3" aperture flux as determined for point 
sources from simulations by Capak et al. (2007). 
To compensate for this approximation, we add in 
quadrature 0.1 and 0.3 mag to the all the measured 
errors in the GALEX NUV and FUV bands. 

For the MIR data we did not degrade the opti- 
cal data to the larger IRAC PSF. Instead, the fol- 



lowing procedures were used. An IRAC flux was 
first associated with the optical sources by match- 
ing their positions in each of the 4 IRAC bands 
within a search radius of 1". Following Surace et 
al. (2004), the IRAC fluxes were then measured 
in a circular aperture of radius 1.9". The aperture 
flux was then converted to a total flux using the 
aperture correction factors 0.76, 0.74, 0.62, 0.58 
at 3.6/im, 4.5/im, 5.6/zm and 8.0/im, respectively 
(Surace ct al. 2004). Since the optical fluxes were 
measured over an aperture of 3" diameter which 
encloses ~ 75% of the flux for a point like source, 
we then multiply the IRAC total fluxes by a fac- 
tor 0.75. This approximation provides good agree- 
ment between the predicted and observed colors 
(z + — 3.6/im and 3. 6/1777 — 4.5/im colors). In order 
to compensate for this approximation, we add in 
quadrature 0.1, 0.1, 0.3 and 0.3 mag to the errors 
in the IRAC 3.6, 4.5, 5.6 and 8.0/im. 

Finally, all magnitudes are corrected for galac- 
tic extinction estimated for each object individ- 
ually, using dust map images from Schlcgcl ct 
al. (1998). We limit the photo-z analysis to 
an area of 2-deg 2 (1.49878 < a < 2.91276 and 
149.41140 < 5 < 150.826934) which has a uni- 
form and deep coverage in all the bands. Poor 
image quality areas (e.g. field boundary, satu- 
rated stars, satellite tracks and image defects) are 
masked. Photo-z are computed only in the non- 
masked regions with a total covered area of 1.73- 
deg 2 . 126 071, 293 627 and 607 617 sources are 
detected at i + < 24, i+ < 25 and i+ < 26. 

2.2. Spectroscopic data 

The spectroscopic samples were observed with 
the VIMOS/VLT spectrograph (zCOSMOS; Lilly 
ct al. 2007) and the DEIMOS/Keck spectrograph 
(Kartaltepe et al. 2008). These two spectroscopic 
samples have very different selection criteria (see 
below); they therefore cover very different ranges 
of redshift and color space, providing a broad sam- 
ple for evaluation of the photo-z. 

The zCOSMOS survey (Lilly et al. 2007) has 
two components: zCOSMOS-bright with a sam- 
ple of 20,000 galaxies selected at i* < 22.5 and 
zCOSMOS-faint with approximately 10,000 galax- 
ies color-selected to lie in the redshift range 1.5 < 
z < 3. In the latter, galaxies are selected by 
color based either on the BzK criterion (Daddi 
ct al. 2004) or the UGR "BM" and "BX" crite- 



4 



filter telescope effective A FWHM s f % at % at % at 

i+ < 24.5 i+ < 25 i+ < 25.5 





CFHT 


3911.0 


538.0 


0.054 


89.3 


85.2 


77.2 


Bj 


Subaru 


4439.6 


806.7 


-0.242 


97.1 


95.2 


90.5 


v 3 


Subaru 


5448.9 


934.8 


-0.094 


99.3 


98.2 


94.2 


g + 


Subaru 


4728.3 


1162.9 


0.024 


96.4 


93.6 


86.0 


r+ 


Subaru 


6231.8 


1348.8 


0.003 


99.6 


99.5 


98.4 


i+ 


Subaru 


7629.1 


1489.4 


0.019 


99.9 


99.9 


99.8 


i* 


CFHT 


7628.9 


1460.0 


-0.007 


37.8 


25.4 


17.4 


z+ 


Subaru 


9021.6 


955.3 


-0.037 


99.8 


97.9 


83.8 


J 


UKIRT 


12444.1 


1558.0 


0.124 


65.4 


49.1 


35.7 


K s 


NOAO 


21434.8 


3115.0 


0.022 


15.3 


10.3 


7.08 


K 


CFHT 


21480.2 


3250.0 


-0.051 


84.1 


68.5 


52.1 


IAA27 


Subaru 


4256.3 


206.5 


0.037 


77.1 


64.3 


48.4 


JA464 


Subaru 


4633.3 


218.0 


0.013 


78.5 


64.3 


47.6 


IAA8A 


Subaru 


4845.9 


228.5 


0.000 


88.7 


78.3 


62.0 


IA505 


Subaru 


5060.7 


230.5 


-0.002 


84.0 


70.0 


52.0 


IA527 


Subaru 


5258.9 


242.0 


0.026 


93.2 


84.6 


68.7 


IA574 


Subaru 


5762.1 


271.5 


0.078 


92.5 


80.0 


60.6 


IAQ2A 


Subaru 


6230.0 


300.5 


0.002 


97.4 


90.2 


72.2 


7A679 


Subaru 


6778.8 


336.0 


-0.181 


99.5 


96.7 


82.9 


IA7QQ 


Subaru 


7070.7 


315.5 


-0.024 


99.8 


97.7 


83.5 


IA738 


Subaru 


7358.7 


323.5 


0.017 


99.5 


94.2 


73.6 


IA7Q7 


Subaru 


7681.2 


364.0 


0.041 


99.6 


93.5 


72.0 


IA827 


Subaru 


8240.9 


343.5 


-0.019 


99.7 


97.2 


81.8 


NB711 


Subaru 


7119.6 


72.5 


0.014 


84.8 


60.7 


41.9 


NB816 


Subaru 


8149.0 


119.5 


0.068 


99.8 


99.1 


88.2 


IRACl 


Spitzer 


35262.5 


7412.0 


0.002 


70.6 


60.8 


48.4 


IRAC2 


Spitzer 


44606.7 


10113.0 


0.000 


62.6 


51.6 


39.7 


IRAC3 


Spitzer 


56764.4 


13499.0 


0.013 


33.7 


26.0 


19.5 


IRACA 


Spitzer 


77030.1 


28397.0 


-0.171 


15.7 


11.3 


8.1 


FUV 


GALEX 


1551.3 


230.8 


0.314 


8.5 


5.8 


4.0 


NUV 


GALEX 


2306.5 


789.1 


-0.022 


19.7 


13.4 


9.2 



Tabic 1: Effective wavelength, width and systematic olfsets sj in magnitude (with our definition, sj have 
the opposite sign to Table. 13 of Capak et al. 2007) plus the fraction of sources in each band with an error 
less than 0.2 magnitude at i + < 24.5, i + < 25, i + < 25.5. 
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rion of Steidel et al. (2004), and the magnitude 
cut was Bj < 24 — 25 (depending on the color 
cut). zCOSMOS-bright galaxies were observed us- 
ing the red grism of VIMOS covering a wavelength 
range 5500A < A < 9000A at a resolution of 600 
(MR grism). For the zCOSMOS-faint sample, ob- 
servations were carried out with the blue grism of 
VIMOS (3600A < A < 6800A) at a resolution of 
200. 

The zCOSMOS-bright survey is now - 50% 
complete. Here we make use of only the extremely 
secure spectro-z measurements with a confidence 
level greater than 99% (class 3 and 4) . This secure 
zCOSMOS-bright sample contains 4148 galaxies 
with a median redshift of ~ 0.48. The zCOSMOS- 
faint survey is in its early stages, and here we use 
a preliminary sample of 148 galaxies with a me- 
dian redshift of z m ~ 2.2 and as faint as i + ~ 25. 
This zCOSMOS-faint spectroscopic sample is not 
fully representative of the average population at 
1.5 < z < 3 due to the selection criteria (e.g. 
Bj < 24- 25). 

The Keck II spectroscopic follow-up of 24^m se- 
lected sources (Kartaltepe et al. 2008) is on-going 
and we refer to this sample as MlPS-spectro-z. 
The DEIMOS spectra cover a wavelength range 
4000A < A < 9000A at a resolution of 600. This 
sample of 24/j,m selected galaxies contains 317 se- 
cure spectro-z (at least two spectral features) with 
an average redshift of z ~ 0.74 and apparent mag- 
nitude in the range 18 < i + < 25. 

For all of the spectroscopic samples used in this 
paper for testing and verification of the photo-z, 
we include only secure spectro-z. Therefore, the 
uncertainties in the spectro-z are neglected and 
the spectro-z are used as a reference to assess the 
quality of the photo-z. 

3. Photometric Redshift derivation 

Photometric redshifts were derived using the 
Le Phare cod^l (S. Arnouts & O. Ilbert) which is 
based on a x 2 template-fitting procedure. In the 
discussion below, we focus on the improvements 
introduced here as compared to Ilbert et al. (2006) 
and the previous COSMOS photo-z (Mobasher et 
al. 2007). 



2 www. oamp.fr/people/arnouts/LE_PHARE. html 



3.1. Galaxy SED template library 




3 3.5 4 4.5 



log ( X(A) ) 

Fig. 1. — SED templates. The flux scale is ar- 
bitrary. The top 12 SEDs (cyan) are generated 
with Bruzual & Chariot (2003). The spiral (green) 
and elliptical (red) SEDs are from Pollctta et al. 
(2007). 

Ilbert et al. (2006) and Mobasher et al. (2007) 
used a set of local galaxy SED templates (CWW: 
Coleman, Wu, & Weedman 1980) which have been 
widely employed for photo-z (e.g. Sawicki 1997; 
Fernandez-Soto et al. 1999; Arnouts et al. 1999; 
Brodwin et al. 2006). Here, we employ a new set 
of templates generated by Polletta et al. (2007) 
with the code GRASIL (Silva et al. 1998). Pol- 
lctta et al. selected their templates for fitting the 
VVDS sources (Le Fevre et al. 2005) from the UV- 
optical (CFHTLS: McCracken et al. 2007) to the 
mid-IR (SWIRE: Lonsdale et al. 2003). There- 
fore, this set of templates provides a better join- 
ing of UV and MIR than those by CWW. The 9 
galaxy templates of Polletta et al. (2007) include 
3 SEDs of elliptical galaxies and 6 templates of 
spiral galaxies (SO, Sa, Sb, Sc, Sd, Sdm). 

We did find that the blue observed colors of the 
spectroscopic sample were not fully reproduced by 
the Polletta et al. (2007) templates. We therefore 
generated 12 additional templates using Bruzual 
& Chariot (2003) models with starburst (SB) ages 
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ranging from 3 Gyr to 0.03 Gyr. We extend the 
BC03 templates beyond 3/ito rest-frame using the 
Sdm template of Polletta et al. (2007). The full 
library of template SEDs, 9 from Polletta ct al. 
(2007) and 12 from Bruzual & Chariot (2003), is 
shown in Fig. [1] Finally, we linearly interpolated 
between some Polletta et al. templates to refine 
the sampling in color — z space. 

Figure [2] shows the observed colors and red- 
shifts of the spectroscopic sample compared with 
the predicted colors for the library SEDs. 

3.2. Emission lines 

Fig. [2] clearly shows how the observed colors os- 
cillate with the redshift, especially when the colors 
are measured with the medium bands (top panels) . 
Comparing the template curves with (dashed) and 
without (solid) emission lines, one sees that the ex- 
pected line fluxes can cause up to 0.4 mag changes 
in the color. This effect is particularly impor- 
tant when colors involving intermediate and nar- 
row band filters are computed (see for example 
the upper-right panel in FigEJ) but can be already 
noticed using broad band colors. The color oscil- 
lations are well explained by the contribution of 
emission lines like Ha, [O III], and [O II] to the 
observed flux; thus the contribution of the emis- 
sion lines to the flux must be taken into account to 
obtain accurate photo-z and this is a major change 
implemented here compared to Ilbert et al. (2006) 
and Mobasher et al. (2007). 

In order to include the emission line contribu- 
tion to the SED, we need to model the emission 
line fluxes (O II, O III, H/3, Ha, Lya) at any 
redshift, template and extinction. The rescaling 
of the template (A in eqnj2]) determines also the 
emission line fluxes (therefore, the modeling of the 
fluxes must be done galaxy per galaxy) . 

Our new procedure estimates the [O II] emis- 
sion line flux from the UV luminosity of the 
rescaled template, using the Kennicutt (1998) 
calibration laws. In the template fitting, a 
UV rest-frame luminosity corrected for dust ex- 
tinction can be computed at every step of the 
redshift /template/extinction grid (the rescal- 
ing factor A is taken into account in the UV 
luminosity). The UV luminosity (at 2300A) 
is then related to the SFR using the relation 
SFR(M Q yr- 1 ) = 1.4 x lO" 28 ^ [erg s^Hz' 1 ) 



from Kennicutt (1998). This SFR can then be 
translated to an [O II] emission line flux using 
the relation SFR (Mq yr^ 1 ) = (1.4 ± 0.4) x 
10" 41 L [OII] (erg s" 1 ) (Kennicutt 1998). This 
translates to : 

log(F [0 u]) = -0.4 x M UV + 10.65 - (1) 

where DM is the distance modulus, Fjon] is ex- 
pressed in units of 10~ 17 erg s _1 cto 2 and M vv 
is the dust corrected Z7U(2300A) absolute magni- 
tude. 




1 2 

-0.4xM(UV) - DM(z)/2.5 + 10.65 



Fig. 3. — Relation between the [O II] flux and the 
absolute mag nitudes in UV (2300A). The solid line 
corresponds to the relation obtained by applying 
the Kennicutt (1998) relations between SFRon 
and SFRuy as used here to include emission 
line fluxes in the photo-z fitting (see eqnQ]) and 
the points are observed emission line fluxes from 
VVDS (Lamareille et al. 2008). (The UV lumi- 
nosities and O II fluxes are corrected for dust ex- 
tinction.) The red and blue points are galaxies 
with Mu-M R > 1.6 and M V -M R < 1.6, respec- 
tively. The larger symbols correspond to larger 
equivalent width. 

Figure [3] shows the measured [O II] fluxes from 
VVDS (Lamareille et al. 2008) and the relation 
(solid line) expected from Kennicutt (1998). We 
can perform this comparison only at 0.5 < z < 1.4, 
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Fig. 2. — The observed colors and redshifts (z s ) for the spectroscopic sample galaxies (open stars). The 
solid lines are the predicted colors as a function of redshift for some SEDs of the library (red: Elliptical and 
green: Spiral from Polletta et al. 2007 ; cyan: Bruzual & Chariot 2003). The solid curves are the predicted 
colors without including emission lines (no reddening for the elliptical templates, E(B — V) = 0.2 for the 
late types) whereas the dashed curves are the same templates including the emission line fluxes (assuming 
Mpuv = —20 in this example). The top right panel clearly shows that the emission lines can change the 
colors up to 0.4 mag. 
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when the [O II] line is observable in the VIMOS 
spectra. Figure [3] shows the good correlation be- 
tween the measured and predicted [O II] fluxes, 
with an rms dispersion of 0.2 dex. In acknowl- 
edgment of this dispersion, we allow the intrinsic 
[O II] flux to vary by a factor of 2 in the tem- 
plates with respect to the nominal flux predicted 
by eqnQ] 

With our procedure, we can predict for each 
galaxy the [O II] flux at every redshift, template 
and extinction combination. For the other emis- 
sion lines, we adopt intrinsic, unextincted flux ra- 
tios of [O III/O II] = 0.36; [H/3/0 II] = 0.61; 
|Ha/0 II] = 1.77 and [Lya/O II] = 2 (McCall 
et al. 1985, Moustakas et al. 2006, Mouhcine et 
al. 2005, Kennicutt 1998). When we apply an 
additional extinction to the template, we modify 
these ratios with the corresponding attenuation. 
Then, we sum the emission line fluxes to the tem- 
plate continuum before integrating through the fil- 
ter transmission curves. 

The effect of these emission lines on the mod- 
eled color — redshift relation is shown in Fig.[2]for 
a galaxy at Mjjv — —20. The oscillations in the 
observed colors versus redshift are well reproduced 
by the models. 

3.3. Systematic offsets 

The x 2 template-fitting method is meaningful 
only if the color — z relation predicted from the 
templates is a good representation of the observed 
color — z relation. Uncertainties in the zero-point 
offsets of photometric bands can lead to systematic 
shifts between the predicted color — z relation and 
the observed colors of the spectroscopic sample. 

To evaluate the zero-point errors, we use the 
spectroscopic sample and set the redshift to the 
spectro-z value. Then, we determine the best-fit 
template for each spectroscopic galaxy. For ran- 
dom, normally distributed uncertainties in the flux 
measurements, AFf should be ~0 (where AFf is 
the average difference between predicted and ob- 
served fluxes in the filter /). Instead, we initially 
find systematic offsets of AFf (as for earlier photo- 
z derivations, e.g. Brodwin et al. 2006; Ilbert et 
al. 2006; Mobasher et al. 2007). Such offsets are 
mainly due to: 

• uncertainties in the absolute calibration of 
the photometric zero-points. 



• uncertainties in the color modeling (filter 
transmission curves, incomplete set of tem- 
plates or an incorrect extinction curve). 

These systematic offsets were removed using 
the iterative procedure detailed in Ilbert et al. 
(2006). For each filter, /, we estimate the values 
s f which minimize AFf . After having applied the 
corrections, s/, in each band /, the systematic off- 
sets derived in a second iteration can change. The 
values converge after 3 iterations and we adopted 
the systematic offsets listed in Table. [TJ 

3.4. Extinction law 

Ilbert et al. (2006) adopted the dust extinction 
law measured in the Small Magellanic Cloud (Pre- 
vot et al. 1984). However, considerable changes in 
the extinction curve are expected from galaxy to 
galaxy. Maraston et al. (2006) considered as a free 
parameter the different extinction curves of the 
"Hyperz" package (Bolzonella et al. 2000) (e.g. 
Milky Way, large and small Magellanic clouds and 
Calzetti). In order to limit the risk of catastrophic 
failures, we adopt an intermediate approach here, 
using the most suitable extinction curve depend- 
ing on the SED template. 

For each galaxy of the zCOSMOS sample, we 
set the redshift to the spectroscopic redshift value. 
Then, we determine the best fit-template and 
the appropriate color excess E(B — V) hest . In 
this fit, we assume an extinction curve k(X) (= 
A(X)/E(B — V)). Since the extinction curves do 
not differ strongly at A > 3000A (blue and green 
curves in Fig. [4]), we fit the templates using only 
passbands with A > 3000(1 + z)k. With this pro- 
cedure, the E(B — V) hcst value does not depend 
significantly on the adopted extinction curve. 

The extinction curves differ strongly at A < 
3000A. Therefore, we use the rest-frame ob- 
served SEDs at Arost-framo < 3000A to discrimi- 
nate between the different extinction curves. Us- 
ing m obs to represent the observed magnitude and 
TO tmcor a * e to represent the predicted magnitude 
from the best-fit template (uncorrected for extinc- 
tion), the extinction A(X) is given by A(X) = 
m ohs - TO*^P lato . Fig. H shows the rest-frame 
(m obs -m t u e n ^P latc )/i;(B-y) bcst for each flux mea- 
surement, i.e the extinction curve k(X). These 
points are compared with the Calzetti et al. (2000) 
and the Prevot et al. (1984) extinction curves. 
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0.2<z<1.3 




3 3.5 4 4.5 3 3.5 4 4.5 



log(X(A)) 



Fig. 4. — Attenuation by dust as a function of A. The points are the extinction Ai/(Ei(B — V)) estimated 
from the galaxies with a spectro-z (see §3.4). The red points at A < 3000A are not used to estimate Ei(B-V). 
The Prevot et al. (1984) and Calzetti et al. (2000) extinction curves are shown with the blue and green solid 
lines, respectively. The Prevot et al. (1984) extinction law is rescaled to the same Ay as the Calzetti law 
by applying a factor 4.05/2.72. The extinction curve derived by Prevot et al. (1984) is used for the galaxies 
redder than SB3 and the Calzetti et al. (2000) extinction law for the galaxies bluer than SB3. 
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(We scaled the Prevot et al. curve to the same 
Ay value as Calzetti et al. by applying a factor 
4.05/2.72 which is the ratio between the Ry of 
the Calzetti and Prevot laws.) The Small Magel- 
lanic Cloud extinction curve (Prevot et al. 1994) is 
well suited for galaxies redder than the starburst 
template SB3. For galaxies bluer than SB3, the 
Calzetti et al. (2000) extinction curve is found 
to be more appropriate (see Fig. |4|). For the IR 
sources which are strongly star-forming, Caputi et 
al. (2008) show also that the Calzetti et al. (2000) 
extinction law is more appropriate. This result is 
not surprising since the Calzetti law was deter- 
mined from observed starburst galaxies. A broad 
absorption excess at 2 175 A (UV bump) seems nec- 
essary to explain the UV flux in some starburst 
galaxies. The presence of this UV bump can be 
seen in a theoretical modeling of the Calzetti law 
(Fischera et al. 2004) and in the K20 sample of 
high redshift galaxies at 1 < z < 2.5 (Noll et al. 
2007). 

To summarize, we apply an additional extinc- 
tion to the templates according to a grid E(B — 
V) = 0,0.05,0.1,0.15,0.2,0.25,0.3,0.4,0.5. We 
use the Prevot et al. (1994) extinction curve for 
the templates redder than SB3, and Calzetti et al. 
(2000) for the templates bluer than SB3. We al- 
low an additional bump at 2175A for the Calzetti 
extinction law if it produces a smaller x 2 ■ No red- 
dening is allowed for galaxies redder than Sb. 

3.5. x 2 minimization 

The photo-z is the redshift value which mini- 
mizes the merit function x 2 (z, T, A): 



x 2-Y / ( F ° hs ~ A xF » Uz,T) 10 ~ 



a 



f 

obs 



(2) 

where F^ Ted (T, z) is the flux predicted for a tem- 
plate T at redshift z. F^ hs is the observed flux and 

^obs ls t ne associated error. The index / refers to 
each specific filter and Sf is the zero-point offset 
listed in Tabled] The opacity of the inter-galactic 
medium (Madau et al. 1995) is taken into account. 
The photo-z is estimated from the minimization 
of x 2 with respect to the free parameters, z, T 
and the normalization factor A. The color excess 
E(B — V) is included in the term T (see section 
3.4). The grid spacing in redshift is Sz = 0.01 and 



the final redshift is derived by parabolic interpola- 
tion of the redshift probability distribution. The 
redshift probability distribution function (PDFz) 
is derived directly from the X 2 ( z ) distribution: 



P(z) cc exp{- 



x 2 W 



x 2 ■ 

A mm 



(3) 



The minimum and maximum redshifts around the 
photo-z solution, corresponding to the la errors, 
are estimated from the equation x 2 { z ) = Xmin + !• 
As in Ilbert et al. (2006), we increased the SEx- 
tractor flux errors by a factor of 1.5. This factor 
does not shift the best-fit photo-z value but broad- 
ens the x 2 peak and derived redshift uncertainty. 



I 2 

+ 




0.5 1 

r + -i + 



Fig. 5. — Color-color diagnostics for sources in- 
cluding stars and galaxies. Stars revealed by the 
X 2 classification are shown in green; the red and 
grey points are galaxies with z p > 1 and z p < 1, 
respectively. 

3.6. Star/AGN/galaxy classification 

For each object, % 2 is evaluated for both 
the galaxy templates and stellar SED templates 
(Chabrier et al. 2000, Bixler et al. 1991). If 

Xgai - Xstar > 0, where 

Xgal an d Xstar arC the 

minimum x 2 values obtained with the galaxy and 
stellar templates respectively, the object is flagged 
as a possible star. Leauthaud et al. (2007) cat- 
alogued point-like sources in the COSMOS field 
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using the peak surface brightness measured on 
the ACS images. The SED \ 2 an d morphological 
classification methods were found to be in excel- 
lent agreement - 84% of the point-like sources at 
i + < 24 are classified as stars with the SED \ 2 
criterion mentioned above, while only 0.2% of the 
extended sources on the ACS images are misclassi- 
fied as stars with the same x 2 criterion. As shown 
in Fig. El the star sequence colors are well distin- 
guished from the bulk of the galaxy population if 
the colors include a NIR band. NIR and MIR data 
are crucial to separate stars from galaxies. There- 
fore, we limit the % 2 star classification to K < 24 
or i^^m > 2% of the point-like sources 

with i + < 24 have K > 24 and F 3 , 6fim > l^Jy. 
These objects, which are only 2% of the total pop- 
ulation at i + < 24, could be stars not recognized 
as such by our SED % 2 classification. 

AGN can be identified by their point-like X- 
ray emission. Most (~ 90%) of the 1887 sources 
detected in XMM-COSMOS (Brusa et al. 2007) 
are dominated by an AGN. Their photo-z deter- 
minations require a different treatment: a correc- 
tion for variability of the photometric data and 
the use of SED templates specifically tuned to 
AGN and their host galaxies. Accurate photo-z 
for the XMM- COSMOS sources are derived in a 
companion paper (Salvato et al. 2008) . Due to the 
flux limit in XMM-COSMOS (Brusa et al. 2007), 
the moderately luminous (log(Lx)=42-43 erg/s) 
and luminous (log(Lx)=43-44 erg/s) AGN are not 
sampled by XMM-COSMOS at z = 0.5 - 1.25 and 
z = 1.5 — 3, respectively. Therefore, this popula- 
tion is not identified as AGN in the galaxy cata- 
logue. 

4. Photometric redshift accuracy 

In this section, we assess the quality of the de- 
rived photometric redshifts by two approaches: 
comparison with high confidence spectroscopic 
redshifts and an analysis of the width of the red- 
shift probability distribution function obtained 
from the \ 2 fitting for the photo-z. The latter 
approach is justified by the excellent agreement of 
the z s -z p distribution with the width of the prob- 
ability distribution function for the spectroscopic 
sample. 



4.1. Comparison of Photometric and Spec- 
troscopic Redshifts 

We first assess the quality of the photo-z by 
comparison with the spectro-z. If Az — z s — z p , we 
can estimate the redshift accuracy from <J& z /{i+z B ) 
using the normalized median absolute deviation 
(NMAD: Hoaglin et al. 1983) defined as 1.48 x 
median ( | z p - z s \ / (1 + z s )) . The NMAD is directly 
comparable to other papers which directly quote 
the rms/(l+z). This dispersion estimate is robust 
with respect to catastrophic errors (i.e., objects 
with \z p — z s \/(l + z s ) > 0.15). The percentage of 
catastrophic errors is denoted by r\. 

Figure [5] (left panel) shows the comparison be- 
tween z p and z s for the zCOSMOS-bright sample 
selected at i + < 22.5. The spectro-z sample is 
selected only by apparent magnitude and is there- 
fore representative of the entire i + < 22.5 pop- 
ulation. We obtain an accuracy of paz/(i+z) = 
0.007 at i + < 22.5 and the distribution of offsets 
((z p — z s )/(\ + z s )) is well fit by a gaussian with 
cr = 0.007 (Fig. [6] right panel). The percentage of 
catastrophic failures is below 1%. 

The zCOSMOS-faint sample (Fig.[0 left panel) 
with a median apparent magnitude of i^ ed ~ 24 
provides a quality check for the photo-z at 1.5 < 
z < 3 where the photo-z are expected to have a 
significantly higher uncertainty. The faint sample 
includes galaxies at i + ~ 25 and the Balmer break 
is shifted into the NIR where the filter set has 
gaps. At 1.5 < z < 3, the accuracy is found to be 
cr Az/(i+z) = 0.06 with 20% catastrophic failures. 
The fraction of high-rcdshift galaxies (z s > 1.5) 
for which a low-rcdshift photo-z (z p < 0.5) was 
assigned is 7% but this failure rate drops to 4% 
if the sample is restricted to galaxies detected 
at IRAC(3.6^m) > lfxJy. The zCOSMOS-faint 
sample actually includes 54 galaxies with low spec- 
troscopic redshift (z s < 0.5 ); 53 out of the 54 
galaxies (> 98%) were assigned the correct photo- 
z low redshift. In summary, we conclude that, for 
the faint galaxies, the photo-z of low redshift ob- 
jects are still assigned correctly, while the failure 
rate for high redshift objects is significantly re- 
duced if the objects have good IRAC detections. 
Such result is expected since the photo-z are al- 
ready including the same information present in 
high-z color selections such as BzK (Daddi et al. 
2004) or similar diagnostics using the IRAC bands 
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z s (z. - z p ) / (l+z s ) 

Fig. 6. — Left panel: Comparison between z p and z s for the bright spectroscopic selected sample 17.5 < 
< 22.5 (zCOSMOS-bright: Lilly et al. 2008). The dotted and dashed lines are for z v = z s ± 0.15(1 + z s ) 
and z p = z B ± 0.05(1 + z s ), respectively. The la dispersion, the fraction of catastrophic failures and the 
median apparent magnitude are listed in the top left corner of the left panel. Right panel: Az/(1 + z s ) 
distribution. The dashed line is a gaussian distribution with a — 0.007. 




Fig. 7. — Comparison between z p and z s . Left panel: zCOSMOS-faint sample (Lilly et al. 2008). The 
open triangles are objects with a secondary peak in the redshift probability distribution function. Right 
panel: infrared selected sample (Kartaltepe et al. 2008) split into a bright sample i + < 22.5 (black), a faint 
sample 22.5 < i + < 24 (red) and a very faint sample 24 < i + < 25 (green). 
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(Fig. [5]) to isolate the good redshift range. 

The zCOSMOS-bright and zCOSMOS-faint 
samples do not probe 0.2 < z < 1.5 at i + > 22.5. 
Here, we use as a comparison sample the MIPS- 
spectro-z (Kartaltepe et al. 2008). In Fig. [7] (right 
panel), we split the MlPS-spectro-z sample into 
bright (i+ < 22.5), faint (22.5 < i+ < 24) and very 
faint samples (24 < i + < 25). For the bright sub- 
sample, the dispersion of the MlPS-selected galax- 
ies is oa 2 /(i+ 2s ) = 0.009, only slightly greater than 
that of the optically selected sample at i + < 22.5. 
For the faint sub-sample (median apparent mag- 
nitude i + ~ 23.1), oaz/(i+z s ) = 0.011, i.e. slightly 
worse than for the brighter optically-selected ob- 
jects. For the very faint sub-sample, the accuracy 
is degraded to aAz/(i+z B ) = 0.053 with a catas- 
trophic failures rate of 20%. This degradation 
is due to decreasing signal-to-noise photometry 
for faint objects and could be amplified by the 
infrared selection which picks up more heavily ob- 
scured galaxies and higher redshift galaxies (e.g. 
Fig.4 of Lc Floc'h et al. 2005). 




1 2 

-z | / (la error) 



Fig. 8. — Cumulative distribution of the ratio \z p — 
z s |/(1(7 error). 65% of the photo-z have a spectro-z 
solution encompassed within the lcr error, close to 
the expected value of 68% (magenta dashed line). 
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Fig. 9. — lcr uncertainty for the z p estimate as a 
function of redshift in different apparent magni- 
tude bin. Each value is computed with 50 galax- 
ies per bin. The thick solid green lines, the solid 
red line, the dashed blue lines and the dotted ma- 



genta lines are for i"" 
24 < i+ < 25, 25 < i< 



< 22.5, 22.5 < i + < 24, 
< 25.5, respectively. 



4.2. Accuracy derived from the Photo-z 
Probability Distribution Function 

Since evaluation of the photo-z accuracy from 
the comparison with spec-z is limited to specific 
ranges of magnitude and redshift, we use the la 
uncertainty in the derived photo-z probability dis- 
tribution to extend the uncertainty estimates over 
the full magnitude/redshift space. 

The reliability of the la uncertainty estimate 
for the photo-z as derived from the Probability 
Distribution Function (PDFz) (see §3.5) can be 
checked by comparing this uncertainty with that 
derived directly from the photo-z - spec-z offsets 
for the spectroscopic sample. Figure [8] shows the 
cumulative distribution of these offsets normal- 
ized by the lcr uncertainty in the probability func- 
tion for the zCOSMOS-bright sample (the ratio 
\z p — z s |/(lcr error) is lower than 1 if the mea- 



sured offset 



is lower than the lcr uncer- 



tainty). 65% of the within the lcr error 

bars, whereas the expected fraction is 68%. We 
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therefore conclude that the la uncertainties in the 
probability function as derived here provide a ro- 
bust assessment of the accuracy in z p . 

Figure [5] shows the la negative and positive un- 
certainties derived from the probability function 
as a function of redshift and apparent magnitude. 
Two clear conclusions emerge: the accuracy is in- 
evitably degraded for fainter galaxies at all red- 
shifts and the photo-z have significantly higher un- 
certainty at z > 1.25. 

From z > out to z = 1.25, the la errors 
do not depend significantly on the redshift, with 
a Az ^ 0.02 at i + < 24 (note that we dropped 
out the division of a by (1 + z) in this anal- 
ysis). To first order, the photo-z are accurate 
when the wavelengths of the Balmer and/or Ly- 
man breaks are well constrained. Therefore, the 
wavelength coverage of the filter set and the pho- 
tometry sensitivities determine the photo-z accu- 
racy as a function of redshift. The COSMOS 
photometry coverage is continuous and dense in 
the optical from the u* band (A e ff ~ 39 11 A) to 
the z + band (A e ff ~ 9021 A). The average wave- 
length spacing between consecutive filters is also 
only 230A. The Balmer break at redshift z = 
out to z = 1.25 is always at A < 9000A, thus the 
high accuracy in this redshift range. 

The accuracy degrades at z > 1.25, where the 
4000A Balmer break goes out of the z + band ; at 
z ~ 1.8 it is a factor 3 larger than at z ~ 1.25 
(caz ~ 0-14 for i + ~ 24). The lack of coverage 
between the z' and the J band accounts for the 
discontinuous increase in uncertainty at z > 1.25 
since the Balmer break can not be located with 
precision. At z > 1.5, the estimated error from 
the photo-z - spectro-z comparison is oAz = 0.19 
(see §4.1; oaz/(i+z s ) = 0-06 at the median redshift 
z ~ 2.2). This accuracy is ~ 1.4 bigger than that 
estimated from the la uncertainty (ctaz ~ 0.14 
at i + ~ 24). Errors due to bias in photo-z will 
not be recovered from the probability distribution 
function which could explain this difference. At 
z > 1.25, all the optical filters, which are about 
80% of all used filters, are sampling the rest-frame 
UV at A < 3500A. Uncertainties in the adopted 
extinction law have considerable impact on the UV 
slope and could introduce small biases in the esti- 
mate of the photo-z. 

At z > 2.5, the accuracy improves improves 
again (ctaz ^ 0.1 for 24 < i + < 25) when the 



4000A Balmer break enters in the J band (z ~ 2) 
and the UV light shortward of L a enters the u* 
band (z ~ 2.3). 

4.3. Importance of zero-point offsets and 
emission lines 

The major improvements in the technique im- 
plemented here are the carefully iterative evalua- 
tion of the photometric zero-point offsets for all 
bands and the allowance for a range of emission 
line contributions to the fluxes of template SEDs. 

The left panel of Fig. ITOl shows the comparison 
between z p and zs computed without correction of 
the systematic band offsets (see §3.3) and without 
including the emission lines (see §3.2). Some sys- 
tematic biases (horizontal and vertical stripes in 
Fig. [T0|) greater than 5 Z ~ 0.1 are clearly seen in 
the photo-z estimate. For example, galaxies with 
z s ~ 0.8 are often shifted to z z ~ 0.7. With the 
iterative calibration of the band offsets turned on 
(right panel of Fig. ITO]) . these biases are limited 
to 5 Z < 0.05. This clearly demonstrates the im- 
portance of zero-point calibration to reduce the 
photo-z bias (as already shown by Brodwin et al. 
2006 and Ilbert et al. 2006). 

The right panel of Fig. [T0l shows the comparison 
between z p and z s without including the emission 
lines (see §3.2). The accuracy is cta 2 /(i+2 s ) ~ 0.02. 
When the emission lines are included in the tem- 
plates, the accuracy is improved by a factor of 2.5 

(0-A*/(i+*e) = °- 007 > left P anel of Fig. E]). There- 
fore, including emission lines in the SEDs is cru- 
cial, especially when the medium bands are used 
to measure z p . 

5. Discussion 

5.1. Comparison with Other Surveys 

We now compare the accuracy of the COSMOS- 
30 photo-z (this work) with those obtained pre- 
viously for COMBO-17 (Wolf et al. 2004) and 
CFHTLS-DEEP (Ilbert et al. 2006). The accu- 
racies can be compared as a function of apparent 
magnitude using the la measured uncertainties in 
the photo-z probabilities. Once again, we check 
also the validity of the la uncertainties for other 
surveys following the method described in §4.2. 
For the COMBO-17 and CFHTLS photo-z, we use 
the spectro-z from the VIMOS-VLT Deep Survey 
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Without band offsets and without emission lines With band offsets and without emission lines 




Fig. 10. — Comparison between z p and z s at different steps of the method. The comparison is done for the 
zCOSMOS-bright sample selected at 17.5 < i + < 22.5. The left panel shows the photo-z computed with a 
standard x 2 method (no emission lines contributions and without calibration of the band offsets). The right 
panel shows z v computed with the calibration of the band offsets, but without including emission lines in 
the templates. Including emission lines improves the accuracy by a factor ~2.5. 



(Le Fevre et al. 2004, 2005). Figure Q3] shows 
that about 53%, 67% of the values for z s are in- 
cluded inside the la uncertainties for COMBO-17 
and CFHTLS-DEEP survey, respectively. The la 
errors for z p in COMBO-17 are rescaled by ~1.2 
to obtain 68% of the z s within the la error (This 
rescale is small and changes nothing in our con- 
clusions). 

Figure [T2] shows the redshift dependence of the 
la uncertainties in photo-z as a function of mag- 
nitude for 0.2 < z < 1.25. The z p for the 
CFHTLS-DEEP were derived from 5 broad bands 
(u*, g', r' , i' , z'). The accuracy of the COS- 
MOS photo-z is improved by a factor of 3 when 
compared to CFHTLS-DEEP. This improvement 
is largely due to the 12 medium bands in this red- 
shift range (NIR data have a real impact only at 
z > 1.25). We checked this by deriving photo-z 
without the 12 medium bands, obtaining an ac- 
curacy of o"a*/(i+zb) ~ 0-03 at i + < 22.5 (similar 
to the COSMOS release of Mobasher et al. 2007 
using 8 broad bands u* , Bj, Vj, g + , r + , i + , z + 
and Ks). COMBO-17 includes 12 medium bands 
in addition to the 5 broad bands and Wolf et al. 



(2004) achieved an accuracy of caz/(i+z) = 0.02 at 
i+ ~ 21.5. The accuracies of COMBO-17 photo- 
z are intermediate between those of COSMOS-30 
and CFHTLS for i\ B < 22.5 but become larger 
that the CFHTLS accuracies at fainter magni- 
tudes, because the COMBO-17 data are about 
1.5 mag shallower than the CFHTLS data. Since 
only secure photo-z (one solution) have been kept 
in the COMBO-17 catalogue, it explains the flat- 
tening of the error bars at i\ B > 24. 

5.2. Redshift Distribution of Galaxies 

Fig. H31 shows the galaxy redshift distribution 
per deg 2 from the COSMOS and CFHTLS sur- 
veys (for 21.5 < i+ uto < 24.5). Although the sur- 
veys are largely independent and the photo-z are 
computed with different codes, the overall agree- 
ment in the redshift distributions is excellent. The 
agreement is particularly good between COSMOS 
and CFHTLS-D2 which covers 1 deg 2 within the 
COSMOS field. 

The density of galaxies at z > 1.5 in the 
CFHTLS fields is of a crucial interest for weak 
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mag range a b c A average z median z 



22.0 < 


i + 

auto 


< 22.5 


0.497±0.019 


12.643±0.409 


0.381±0.016 


4068.19 


0.66 


0.66 


22.5 < 


7 + 
auto 


< 23.0 


0.448±0.016 


9.251±0.218 


0.742±0.030 


9151.98 


0.76 


0.72 


23.0 < 


?+ 

auto 


< 23.5 


0.372±0.012 


6.736±0.094 


1.392±0.055 


18232.24 


0.90 


0.82 


23.5 < 


7 + 
auto 


< 24.0 


0.273±0.008 


5.281±0.039 


2.614±0.096 


35508.58 


1.05 


0.92 


24.0 < 


7 + 
auto 


< 24.5 


0.201±0.005 


4.494±0.024 


3.932±0.134 


60306.30 


1.18 


1.00 


24.5 < 


7 + 
auto 


< 25.0 


0.126±0.003 


4.146±0.021 


5.925±0.191 


103340.04 


1.25 


1.06 



Table 2: Galaxy redshift distribution per deg 2 . The parameters a, b, c and A of the cqn|4] function (Fu et 
al. 2008) are given per apparent bin. 




Fig. 11.— Same as Fig[5]for the CFHTLS-Deep 
survey (Ilbert et al. 2006) and the COMBO-17 
survey (Wolf et al. 2004). 

lensing analysis (Benjamin et al. 2007). The 
CFHTLS redshift distribution presents a bump at 
z ~ 3, but this excess could be due to misidentifi- 
cations between the z < 0.4 and z > 1.5 photo-z 
(Van Waerbeke et al. 2008) when no NIR data are 
available (as is the case for CFHTLS). The COS- 
MOS photo-z are computed with NIR data and 
such catastrophic failures are limited (see §4.1). 
Fig. Q2] shows that this bump seems created by a 
deficit of galaxies at 2.5 < z < 3, by comparing 
CFHTLS-D2 and COSMOS. 



Fig. 12. — la error for the z p estimate as a func- 
tion of the apparent magnitude in the redshift 
range 0.2 < z < 1.25. The la errors have been 
rescaled by a factor 1.2 for the COMBO-17 survey 
(Wolf et al. 2004) (see text and Fig.fTTj). 

Fu et al. (2008) used the CFHTLS photo-z dis- 
tribution to estimate the matter density param- 
eter fi m and the amplitude of the matter power 
spectrum a&. The redshift distribution from Fu et 
al. (2008) at 21.5 < i+ < 24.5 and z < 2.5 is over- 
plotted in Fig. [Rj] with and without (J. Coupon, 
private communication) the weight applied for the 
weak lensing selection (dashed and dashed-dot 
lines, respectively). Given the 20% variation ex- 
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Fig. 13. — Redshift distribution from the 2-deg 2 
COSMOS survey (black solid line), from the 1- 
deg 2 CFHTLS-D2 field (green dashed line), from 
the 4-deg 2 CFHTLS-DEEP fields (blue dotted 
line, including also D2). CFHTLS-D2 covers 1 
deg 2 within the COSMOS field. The red long 
dashed line is the redshift distribution obtained 
by Fu et al. (2008) who fit the CFHTLS-DEEP 
photo-z in the magnitude bin 21.5 < i + < 24.5 
and the dashed-dotted line is obtained without the 
weight applied for the CFHTLS weak lensing se- 
lection (J. Coupon, private communication). 

pected from cosmic variance (Ilbert et al. 2006), 
it is in excellent agreement with the COSMOS- 
30 redshift distribution. This agreement suggests 
that the derivation of cr 8 (f7 m /0.25) ' 64 = 0.785 ± 
0.043 by Fu et al. (2008) is not suffering from bi- 
ases due to the photo-z or significant cosmic vari- 
ance. 

The COSMOS-30 redshift distribution can be 
fit with a parametrization similar to that used by 
Fu et al. (2008): 



/ x Az a + Z ab ) 

(z° + c) 



(4) 



where A is the normalization factor and a, b and 
c are free parameters. The best fit parameters 
a, b and c are given in Table [21 as well as the 



Fig. 14. — Top panel: Evolution of the redshift 
distribution as a function of i + magnitude in the 
COSMOS field. We use the parametrization of Fu 
et al. (2008) in different redshift bins. Bottom 
panel: Cumulative redshift distribution. 

median redshifts. Fig. Q3] shows the best fit red- 
shift distribution per apparent magnitude bin. As 
expected, the median redshift increases at fainter 
apparent magnitude, ranging from z m = 0.66 at 



22 < i+ < 22.5 to 



6. Summary 



1.06 at 24.5 < i+ < 25. 



This paper presents a new version of the pho- 
tometric redshift catalog for the 2-deg 2 COS- 
MOS survey computed with new ground-based 
NIR data, deeper IRAC data and a new set of 
12 medium bands from the Subaru Telescope. 
The COSMOS photometry now includes a to- 
tal of 30 filters - from the UV ( GALEX) to the 
MIR (Spitzer-IRAC) . The photo-z catalogue de- 
rived here contains 607,617 sources at i + < 26. 
The 1887 XMM-COSMOS sources (mainly AGN) 
are not included in this catalogue; their photo-z 
are derived in Salvato et al. (2008) with simi- 
larly good accuracy using a set of templates for 
composite AGN/galaxies. 

The galaxy photo-z were tested and improved 
using spectroscopic redshift samples from the 
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Sample 



N median z s median i + &Az/(i+z s ) V m % 



zCOSMOS bright 


17.5 < i 


f < 22.5 


4146 


0.48 


21.6 


0.007 


0.7 


zCOSMOS faint 


1.5 < z s 


< 3 


147 


2.20 


24.0 


0.054 


20.4 


MIPS bright 


17.5 < i 


f < 22.5 


186 


0.68 


21.7 


0.009 


0.0 


MIPS faint 


22.5 < i 


f < 24.0 


116 


0.90 


23.1 


0.011 


0.0 


MIPS very faint 


24.0 < i 


f < 25.0 


15 


1.15 


24.2 


0.053 


20.0 



Table 3: Redshift accuracies estimated from the comparison between photometric redshifts and spectroscopic 
redshifts. 



zCOSMOS and Keck surveys. Biases in the photo- 
z were removed by iterative calibration of the pho- 
tometric band zero-points. As suggested by the 
data, two different dust extinction laws were ap- 
plied specific to the different SED templates. A 
new method to take into account the emission lines 
was implemented using relations between the UV 
continuum and the emission line fluxes associated 
with star formation activity. The allowance for 
emission lines decreased the photo-z dispersion by 
a factor 2.5. 

Based on a comparison between our new values 
for z p and 4148 measured values of z s from zCOS- 
MOS, we estimate an accuracy of <Jaz/(i+z s ) — 
0.007 for the galaxies brighter than i + = 22.5. 
The accuracies measured with the various spec- 
troscopic samples are summarized in table [3] We 
extrapolate this result to fainter magnitudes us- 
ing the la uncertainties in the photometric red- 
shift probability functions. This is found to pro- 
vide reliable uncertainty estimates for the photo-z 
technique developed here since these uncertain- 
ties agree with the dispersion in the offsets of 
photo-z from spec-z for the spectroscopic sample. 
At z < 1.25, we estimate a photo-z accuracy of 
a Az = 0.02, 0.04, 0.07 for i+ ~ 24, i+ ~ 25, 
i + ~ 25.5. The accuracy is strongly degraded at 
i+ > 25.5 and the exploitation of the COSMOS- 
30 photo-z at fainter magnitudes should be done 
carefully. The accuracy is 3-5 times better than 
the photo-z determined for the CFHTLS-DEEP 
(Ilbert et al. 2006) and the previous COSMOS 
photo-z release (Mobasher et al. 2007), and 2 
times better than the accuracy of the photo-z de- 
termined for COMBO- 17 at i+ < 23 (Wolf et al. 
2004). Deep NIR (J, K) and IRAC data were 
essential to keep the catastrophic failures low at 



z > 1.25. We note that the accuracy of the COS- 
MOS photo-z could soon be further improved at 
z > 1.25 with the addition of new data currently 
being obtained in the Y band (UKIRT), H band 
(CFHT) and the ULTRA- VISTA survey. 

Our photo-z catalogue contains 607,617 sources 
at i + < 26. The accurate photo-z derived here for 
this extremely large sample of galaxies are cru- 
cial to scientifically exploit the full legacy value of 
the multi-A data sets in the COSMOS field (HST- 
ACS, Spitzer, GALEX, VLA, XMM&nd Chandra). 
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HST-GO-09822. French co-authors acknowledge 
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project). More information on the COSMOS sur- 
vey is available at http:/ /www. astro. caltech.edu/cosmos 
. We also greatly appreciate the hospitality pro- 
vided by the Aspen Center for Physics where this 
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